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We study the renormalization group flows of the two terminal conductance of a superconducting 
junction of two Luttinger liquid wires. We compute the power laws associated with the renormal- 
ization group flow around the various fixed points of this system using the generators of the SU (4) 
group to generate the appropriate parameterization of an S-matrix representing small deviations 
from a given fixed point §-matrix (obtained earlier in Phys. Rev. B 77, 155418 (2008)), and we 
then perform a comprehensive stability analysis. In particular, for the non-trivial fixed point which 
has intermediate values of transmission, reflection, Andreev reflection and crossed Andreev reflec- 
tion, we show that there are eleven independent directions in which the system can be perturbed, 
which are relevant or irrelevant, and five directions which are marginal. We obtain power laws 
associated with these relevant and irrelevant perturbations. Unlike the case of the two-wire charge- 
conserving junction, here we show that there are power laws which are non-linear functions of V(0) 
and V(2kp) (where V(k) represents the Fourier transform of the inter-electron interaction potential 
at momentum k). We also obtain the power law dependence of linear response conductance on 
voltage bias or temperature around this fixed point. 
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Electron-electron (e-e) interactions in low-dimensional 
systems (one-dimensional (1-D) quantum wires (QW) 
and dots) can lead to non-trivial low energy transport 
properties due to the Luttinger liquid (LL) ground state 
of the system. In this context, a geometry which has 
gained considerable attention in the recent past is the 
multiple LL wire junction. In general, junctions of multi- 
ple QW can be viewed as quantum impurities in a LL from 
which electrons get scattered at the junction. For the 
simplest case of two-wires, the junction can be modeled 
as a back-scatterer while for the general case of multiple 
QW, the junction represents a more non-trivial quantum 
impurity which may not be as straightforward to model 
microscopically. 

For the two-wire system, it is well-known^ 2 - that in 
the presence of a scatterer, there are only two low en- 
ergy fixed points - (i) the disconnected fixed point with 
no transmission (i.e. the transmission amplitude for in- 
cident electron or hole, t — 0) which is stable and (ii) 
the transmitting fixed point with no reflection (t = 1) 
which is unstable. More recently, the low energy dynam- 
ics of multiple LL wires connected to a junction have also 
been studied in detail3^£&Zi£&^ and several interesting 
fixed points have been found, including continuous one- 
parameter families of fixed points 11 . These studies have 
also been generalized theoreticallyi^^^^^^^^ 
to describe a junction of 1-D wires with superconduc- 
tors and have also been generalized to include spin 21 . 
Our recent wor k 22 : 23 ' 24 has generalized these studies to 
the case of superconducting junction of multiple QW. In 
such a system, due to the proximity of the superconduc- 
tor, both electrons as well as holes take part in the trans- 
port which leads to very interesting transport properties 
at small bias resulting from the interplay of LL correla- 



tions and the proximity induced pair potential. 

In this article, we study the case where two LL QW are 
coupled simultaneously to a bulk superconductor. The 
physical separation between the junctions of the two- 
wires with the superconductor is of the order of the 
size of the Cooper-pair. This leads to the realiza- 
tion of a normal-superconductor-normal (NSN) junc- 
tion which allows for direct tunneling of electrons from 
one wire to the other and also allows a finite am- 
plitude for the crossed Andreev reflection (CAR) pro- 
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reflection and Andreev reflection (AR) processes. 

In an earlier study of the NSN junction^ 3 -, we showed 
that the NSN junction has more than two fixed points un- 
like the normal two-wire junction (as mentioned above) 
or thc junction of LL with a bulk superconductor 
(NS junction) which has only two fixed points - (i) the 
Andreev fixed point where the amplitude for Andreev 
reflection (AR), r a = 1 and normal reflection amplitude, 
r = and which is unstable and (ii) the disconnected 
fixed point where ta = and r = 1, and which is 
stabl o 12 : 13 . We showed that there exists a fixed point 
with intermediate values of transmission and reflection. 
Thus, the NSN junction is the minimum configuration 
which possesses non-trivial fixed points with intermediate 
transmission and reflection amplitudes. In what follows, 
we will focus mainly on the NSN junction. 

In the previous studies, a comprehensive analysis of 
the various possible perturbations allowed by symmetry 
around all the fixed points of the NSN junction was lack- 
ing. In this article, we carry out a systematic stability 
analysis for each of the fixed points obtained earlier in 
Ref. [23| for the NSN junction and we predict the power 
laws associated with all possible independent perturba- 
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tions that can be switched on around these fixed points. 
Our analysis provides us with renormalized values of the 
various transmission and reflection amplitudes around 
these fixed point values which can then be used to obtain 
the Landauer-Biittiker conductances. 

We start with a brief review of the RG method followed 
in Refs. [U and where an S- matrix formulation was used 
to compute the linear conductance and inter-electron in- 
teractions inside the QW were taken into account by al- 
lowing the S-matrix to flow as a function of the relevant 
energy scale (like temperature, bias voltage or system 
size) using an RG procedure. This method works well 
when e-e interaction strength inside the QW is weak so 
that it can be treated perturbatively. This is usually re- 
ferred to as the weak interaction renormalization group 
(WIRG) procedure^. 

To benchmark our calculation with known results, we 
first calculate the complete set of all possible power laws 
associated with the independent perturbations that can 
be switched on around the time-reversal symmetry bro- 
ken chiral fixed points (CFP) and the time reversal sym- 
metric Griffiths fixed point (GFP) of a normal junction 
of three LL wiresi Then we apply the same procedure 
to the case of the NSN junction. The strength of our 
formulation to obtain the power law scaling of the per- 
turbations turned on around the various fixed points lies 
in the fact that the same method is applicable to both 
normal as well as superconducting junctions of any num- 
ber of QW. 

Now, for the case of three LL wires meeting at a nor- 
mal junction, let us assume that the wires are parame- 
terized by spatial coordinates Xi which go from zero to 
infinity with the junction being situated at Xi — (i be- 
ing the wire index). The junction can be parameterized 
by a 3 x 3 S-matrix with diagonal elements ru and off- 
diagonal elements Uj. In the presence of e-e interactions, 
the RG equations can be derived 2 ^ by first expanding the 
electron wave- function in each of the wires in terms of re- 
flected and transmitted electron waves (scattering wave 
basis) . Then the amplitude of scattering of electrons from 
the Friedel oscillations in the wires can be deduced by us- 
ing a Hartree-Fock decomposition of the e-e interaction 
term. Finally the RG equation is obtained by applying 
the poor-man's scaling approach^. The RG equations 
for the entire S-matrix can be written in a concise and 
compact from given by 4 
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where, I = \n{L/d) is the dimension- less RG scale (L cor- 
responds to the physical length scale or energy scale at 
which we are probing the system and d is the short dis- 
tance or high energy cut-off). Here F is a diagonal matrix 
with Fa = —aru/2 and a is the repulsive e-e interaction 
parameter which is related to the LL parameter K as 
K = ((l-a)/(l + a)) 1/2 . 

Analogously, the NSN junction can be described in 
terms of an S-matrix with elements describing trans- 
mission of both electrons and holes and their mixing at 



the junction. The corresponding RG equation for the S- 
matrix was obtained by the present author a 22 i 23 which 
was an extension of Eq. [T] to the superconducting case. 
Here too for the NSN case, we will assume that the two 
wires are parameterized by spatial coordinates x\ and 
X2 where xi,X2 vary from zero to infinity and that the 
junction is at X\ — = x%. The presence of the su- 
perconductor is encoded in the parametrization of the 
S-matrix representing the junction. Of course, this way 
of accounting for the presence of the superconductor as- 
sumes that the superconductor imposes static boundary 
conditions on the two wires forming the NSN junction. 
This is a valid approximation as long as one is focusing on 
sub-gap transport properties of the junction. We also as- 
sume that the superconductor at the junction is a singlet 
superconductor; hence the spin of the incident electron 
or hole is conserved as it scatters off the junction. This 
results in a block diagonal form of the S-matrix with each 
spin block being a 4 x 4 matrix representing scattering 
of electrons and holes within the given spin sector. 

The S-matrix at the superconducting junction for the 
spin-up, spin-down, electron-hole and left-right symmet- 
ric (symmetry in wire index) case (suppressing the wire 
index) can be parameterized by r, the normal reflec- 
tion amplitude, r A , the AR amplitude, t A , the CAR am- 
plitude 2 ^^ and t, the transmission amplitude for both 
electrons and holes. The fermion fields can then be ex- 
panded around left and right Fermi points on each wire 
as ip ls (x) = i S (x)e lkpX + ^ is (x) e~ l kp x ; where i is 
the wire index, s is the spin index which can be "f , j and 
1(0) stands for incoming (outgoing) fields. Note that 
^i(o)i x ) are slowly varying fields on the scale of kp 1 . 
Electrons with momenta k in vicinity of hp, on each wire 
at position x is given by 



dk 
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where b ks is the particle destruction operator and d ks is 
the hole destruction operator and we have allowed for 
non-conservation of charge due to the proximity effect 
induced by the superconductor. We then allow for short- 
range density-density interactions between the fermions, 
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Pis') , (3) 



Following the procedure outlined in Ref. [23|, we find that 
the RG equation for the S-matrix continues to be of the 
form given in Eq. [TJ but now F is a non-diagonal matrix, 



(4) 



where a and a' are the cc-independent part of the mean 
field amplitudes for Friedel oscillations and the proximity 
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induced pair potential inside the QW respectively. Gen- 
eralization to particle-hole non-symmetric situations will 
make this matrix asymmetric. It is worth pointing out 
that even though the expression for electron field in Eq.[2] 
assumes particle-hole symmetry which leads to consider- 
able simplification in the derivation for the RG equation 
(Eq. [!}, our formalism is more general. The RG equation 
(with appropriate modification of the F matrix for the 
asymmetric case) will also hold for S-matrices represent- 
ing situations where the wire index symmetry as well as 
the particle-hole symmetry is broken. 

We will mainly focus on three different fixed points 
- the CFP and GFP of a normal junction of three 
LL wires^l and the symmetric fixed point (SFP)22 of 
the NSN junction. First we discuss the stability around 
the CFP and GFP of a normal junction of three LL wires 
(Y-junction) to benchmark our calculation with known 
results^. As a first step towards performing a systematic 
stability analysis, we need to obtain an S-matrix which 
results from a very small unitary deviation from the fixed 
point S-matrix. Given the number of independent pa- 
rameters of the S-matrix dictated by symmetry and uni- 
tarity constraints, the most general deviation from the 
fixed point S-matrix can be obtained by multiplying the 
fixed point S-matrix by another unitary matrix which 
is such that it allows for a straightforward expansion in 
terms of small parameters around the identity matrix. 
This is realized as follows - 



where F is the diagonal part of the following quantity 
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where So represents the fixed point S-matrix and Aj's 
(along with the identity Ao = /) are the eight generators 
of the SU(3) group which are traceless hermitian ma- 
trices. This can be straightforwardly generalized to the 
case of N wires by using SU(N) matrices. Perturbations 
around these fixed points are characterized in terms of 
the ej's. Of course, the resulting S-matrix obtained in 
this way corresponds to a small unitary deviation when 
£j's are small. To first order in fj's, Eq. [5] reduces to 



So [l + i^Xj 
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where So for CFP and GFP fixed points are 
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respectively. Using Eq. the RG equation (Eq. [1} with 
S expanded to the linear order in ej becomes 
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F = 
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Restricting the RHS of Eq. [5] to linear order in ey's, one 
then obtains nine coupled linear differential equations. 
Next by applying a unitary rotation, we can decouple 
these coupled equations (Eq.[5]) and re-cast them in terms 
of new variables e^- (which are linear combinations of the 
original e j ) . The equations are now given by 
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where fij is a real number corresponding to the 'power 
law' associated with perturbations turned on along each 
of the new nine eigen-directions e'j . fij < indicates that 
the given direction is stable and fij > indicates that it 
is unstable. Here the non-diagonal Cj are related to the 
diagonal e'j by ej = J2i ^ji e 'i where U is the diagonalizing 
rotation matrix. 

Hence we obtain all the power laws associated with 
the independent perturbations that can be switched on 
around a given fixed point S-matrix. Now it is straight- 
forward to show that the power laws associated with the 
CFP and GFP are given by [a/2, a/2, 0, a/2, a/2, a/2, 
a/2, 0, 0] and [0, 0, 0, 0, 0, -a/3, 2a/3, 2a/3, a] re- 
spectively which is consistent with results obtained in 
Ref. |^.. The value zero corresponds to marginal di- 
rections while the values with positive or negative signs 
correspond to stable or unstable directions respectively. 
We do not write the explicit form of the U matrix for 
the CFP and GFP as they are needed only for obtaining 
the explicit form of the power law correction to the fixed 
point conductance which we do not calculate for these 
cases. 

Finally, let us discuss the stability around the different 
RG fixed points of the NSN junction. First we focus on 
the SFP of the NSN junction. In the presence of e-e in- 
teraction inside the QW , the incident electron (hole) not 
only scatters from the Friedel oscillations as an electron 
(hole) but also scatters from the proximity induced pair 
potential inside the QW as a hole (electron). Now the 
amplitude of both of these scattering processes are pro- 
portional to the e-e interaction strength inside the QW . 
The competition between these two scattering processes 
which actually arise due the same e-e interaction strength 
inside the QW leads to the presence of the new SFP where 
all the scattering amplitudes have intermediate non-zero 
values. This fact is unique about this fixed point and 
hence this fixed point is the central focus of our discus- 
sion. Details of this fixed point arc further elaborated in 
the discussion at the end of this article. 

We adopt the same procedure as described above for 
the three-wire junction but now with SU (4) generators. 
This is so because the full 8x8 S-matrix describing the 
NSN junction has a block diagonal form with each spin 
block (up and down spin sectors) being represented by a 
4x4 matrix. Hence we have a unitary starting S-matrix 
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deviating from the fixed point S-matrix (So), as given 
before by Eq.[Sl except that now the sum over j runs from 
1 to 16 since Ay's now represent the fifteen generators of 
the ST/ (4) group along with the identity matrix. The So 
which describes the SFP— is given by r = 1/2, t = 1/2, 
v a = —1/2 and tjy = 1/2. Note that the SFP is a particle- 
hole, left-right symmetric fixed point and hence the entire 
4x4 S-matrix is determined completely by the above 
given four amplitudes for r, t, r a^a- 

We then solve Eq. [5] for this case with sixteen cou- 
pled equations up to the first order in the small pertur- 
bations ej's. We obtain the sixteen eigenvalues which 
correspond to the power laws around the sixteen eigen- 
directions. These power laws around the various eigen- 
directions can be listed as [0, 0, 0, 0, 0, —a/2, —a/2, 
(a - a')/2, a'/2, a'/2, (-a + a')/2, (-a + a' )/2, 
(a + a')/2, (a + a') /2, (a - a' - V9a 2 + 14aa / + 9a' 2 ) /4, 
(a - a' + V9a 2 + 14aa' + 9a'")/4]. 
Hence we note that there are five marginal directions, 
two stable directions, four unstable directions and four 
other directions whose stability depends on the sign of 
a — a'. One of the most striking outcomes of this anal- 
ysis is the fact that we obtain two power laws which are 
not just simple linear combinations of V(0) and V(2kp). 
Instead, they appear as square roots of quadratic sum of 
these quantities. Our analysis actually leads to the first 
demonstration of the existence of such power laws in the 
context of quantum impurity problems in LL theory and 
this is the central result of this article. 

Having obtained the power laws the next task is to ob- 
tain an explicit expression for the Landauer-Biittiker con- 
ductance corresponding to perturbations around these 



fixed points along some of the eigen-directions. Now note 
that the RG equation is expressed in terms of e"s whereas 
the S-matrix representing small deviations from the fixed 
point is expressed in terms of e's. The two terminal lin- 
ear conductance across the junction depends explicitly 
on the S-matrix element which are expressed in terms 
of e's (see Eq. Hence in order to obtain an expres- 
sion for conductance in terms of the temperature or the 
applied voltage dependence induced by e-e interaction, 
we need to first assign bare values to the various pertur- 
bations parameterized by e's and then express the e"s 
evolved under RG flow in terms of these bare values of 
e"s as e'(A) = (A/A ) A1 e where A corresponds to the 
energy scale at which we are probing the system (which 
can be either voltage bias at zero temperature or temper- 
ature at vanishing bias) and Ao is the high energy cut- 
off expressed in terms of voltage or temperature. Then 
by using the rotation matrix which diagonalizes the cou- 
pled RG equations, we express e's in terms of e"s written 
explicitly as a function of temperature or voltage. Fi- 
nally plugging these renormalized values of e's into the 
S-matrix given by Eq. O we get all the transmission and 
reflection amplitudes for the system as explicit functions 
of the temperature or voltage carrying the specific power 
laws associated with perturbations switched on along the 
eigen-directions. These amplitudes are now directly re- 
lated to the linear conductances. 

Now we will calculate expression for conductance for a 
simple case where only one of the e'(= e' 15 ) is turned on. 
For this we need the U matrix for this case which is given 
by 
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where, A = 4(a + a')/(a - a' - V9a 2 + 14aa' + 9a' 2 )/4 
and B — 4(a + a') /(a — a' + V9a 2 + 14aa' + 9a'") /4. 
We choose this specific direction to perturb the system 



as this corresponds to a power law which is not a linear 
function of V(0) and V(2Kp) and hence interesting to 
study. The S-matrix to quadratic order in e' 15 is given by 
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So, the scaling of sub-gap conductance (to 0(ef 5 )) for 
an incident electron and a hole taking into account both 
spin-up and spin-down contributions in units of 2e 2 /h is 
given by 
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where e' 15 = e; 5i0 (A/A )( Q - Q '-^ 9aa+14aa ' +9Q ' 2 )/ 4 . Here 
the superscripts e and h stand for electrons and holes 
while the subscripts 1 and 2 stand for first and second 
wire respectively. Also, Gf 2 = \t e A h 12 \ 2 — l*!l| 2 where t ee 
is the transmission amplitude for electrons and t e A rep- 
resents CAR amplitude for electrons. Similar expressions 
hold for the holes. In the expressions of power laws given 
above, a = (32 — 2gi) /2ttHvf and a' = (gi + g-^j^-nhvp 
where the bare values of gi(d) — V(2kp) and g2(d) = 
V"(0). In our stability analysis, we have assumed a < a' 
which is consistent with experimental observations^. For 
the special case when gi = 2<?i, a vanishes and only a 1 
survives. 

It is very interesting to note that even though the S- 
matrix corresponding to perturbation along e' 15 breaks 
both time reversal and electron-hole symmetry, the two 
terminal linear conductance restores particle-hole sym- 
metry. Secondly it might be of interest to note the fact 
that the fixed point conductance admits correction along 
e' 15 which is linear in e' 15 and not quadratic. Normally 
when we perform a stability analysis around a fixed point 
S- matrix whose elements are constituted out of unimodu- 
lar numbers (representing disconnected or perfectly con- 
nected fixed points), it is always possible to identify var- 
ious terms of the §-matrix, representing small unitary 
deviations from the fixed point §o _ma trix in terms of var- 
ious tunneling operators which are perturbatively turned 
on around the fixed point Hamiltonian. Hence a straight 
forward perturbative linear conductance calculation us- 
ing the Hamiltonian along with the tunneling parts will 
suggest that the correction due to the S-matrix repre- 
senting small deviation from fixed point So-matrix must 
introduce correction to fixed point conductance which are 
quadratic in terms of the deviation parameter. But this 



argument applies only to those fixed points which cor- 
respond to completely connected or disconnected wires 
and not to fixed points which have intermediate values 
for various transmission and reflection amplitudes like 
the SFP . In other words, an arbitrary deviation from 
SFP may not be easily representable as a tunneling op- 
erator. This explains why the linear dependence of the 
conductance on e' and hence the corresponding power 
law dependence looks unconventional. 

As a cross check, we see that we get back the power 
laws associated with the symmetric fixed points of the 
four-wire junction once we substitute a' = in the ex- 
pression for the power laws of the SFP for the NSN junc- 
tion. Although our geometry does not correspond to the 
real junction of four LL wires, the presence of both elec- 
tron and hole channel mimics the situation of a four-wire 
junction. More specifically, the symmetric fixed point of 
the NSN junction (SFP) turns out to be identical to the 
symmetric fixed point of the four-wire junction due to 
perfect particle-hole symmetry of the SFP when a' is set 
to zero. 

Next we enumerate and discuss the stability of the 
other fixed points (RFP, AFP, TFP and CAFP) obtained 
in Ref. [HI for the NSN junction : 

(a) t = t A = r A = 0,r = 1 (RFP) : This fixed 
point turns out to be stable against perturbations 
in all directions. There are ten directions for which 
the exponent is -a while two others with the expo- 
nents -(a+a'). The remaining four directions are 
marginal. 

(b) t = t A = r = 0, r A = 1 (AFP) : This is unstable 
against perturbations in twelve directions. There 
are ten directions with exponent a and two direc- 
tions with exponent (a+a'). The remaining four 
directions are marginal, as for RFP. 

(c) r A = t A = r = 0,t = 1 (TFP) : This fixed 
point has four unstable directions with exponent 
a, two stable directions with the exponent -a' and 
the remaining directions are marginal. 

(d) r A = t = r = 0, t A = 1 (CAFP) : This has four 
unstable directions with exponent a and two stable 
directions with the exponent -a' and the remaining 
directions are marginal. 
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Note the close similarity in stability between CAFP and 
TFP fixed points. This can be attributed to the fact 
that both these fixed points belong to the continuous 
family of marginal fixed points defined by the condition 
\t\ 2 + \tA\ 2 = 1. The entire family of fixed points is 
marginal because for these fixed points, the amplitudes 
for Friedel oscillation and pair potential in the wire van- 
ish identically. 

Hence, we notice that for the AFP only the scattering 
amplitude from the pair potential inside the QW is non- 
zero as only r& is nonzero, and for RFP only the scat- 
tering amplitude from Friedel oscillations are non-zero 
as only r is nonzero. Furthermore, both for CAFP and 
TFP , the amplitude for scattering from the Friedel oscil- 
lations as well as from the pair potential is zero as in these 
cases both r and ta are zero. So SFP is the only fixed 
point for which both the amplitude for scattering from 
the Friedel oscillations and the pair amplitude are finite; 
hence, this fixed point is nontrivial. Its very existence 
can be attributed to the interplay of these two differ- 
ent scattering processes arising from Friedel oscillations 
and the pair potential inside the wire. The conductance 
at this fixed point gets contribution from both the elas- 
tic co-tunneling (CT) of electrons through the supercon- 
ductor as well as through the crossed Andreev reflection 
(CAR) process. Since both electron and hole channels 
contribute with opposite signs to conductance, if we give 
a small perturbation around this fixed point, we get an 
interesting non-monotonic behavior of the conductance 
Gnsn = Gcar — Gct • This effect emerges due to the 
competition between the electron and the hole channel 
and it can be of interest from an experimental point of 
view. Also note that at the SFP , CT amplitude of elec- 
trons t = 1/2 and the CAR amplitude t A = 1/2. This 
means that if we have an incident spin-polarized beam of 
(say "up" polarized) electrons on the junction, when the 
junction is tuned to this fixed point, 25% of the spin up 



electrons get transmitted through the junction and 25% 
of the spin up electrons get converted to spin up holes 
as they pass through the junction. Hence the transmit- 
ted charge across the junction is zero on the average, but 
there is pure spin current flowing out of the junction. 
Equivalently, we can think that the pure spin current is 
generated due to flow of two beams of electrons of equal 
intensity, one with spin up electrons and the other with 
spin down electrons propagating in opposite directions. 
Therefore the SFP can be relevant for future spintronics 
applications. These points have been discussed in detail 
in Ref. [H 

To summarize, we have laid down a scheme to per- 
form a systematic stability analysis which works well for 
both normal and superconducting junctions of multiple 
LL QW. Using our procedure, we reproduced the known 
power laws for the fixed points of the three- and four- wire 
junctions. Then we applied it to the NSN junction and 
established the existence of non-trivial power laws which 
are non-linear functions of V(0) and V(2kp)- Finally, 
we calculated the Landauer-Biittiker conductance asso- 
ciated with the perturbations switched on around these 
fixed points and found the explicit voltage or tempera- 
ture power law dependence. 
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There is a correction to the power law obtained for the 
GFP for the three-wire junction in Ref. Q. Lai et. al. had 
predicted a stable direction with power law —a which 
should be corrected to — a/3 3 - 9 - as is obtained in this work. 
The corrected power law is also consistent with that ob- 
tained in Ref. [4TJ using a functional RG procedure. 



